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An intriguing possibility for TeV scale physics is the existence of neutral long lived particles 
(LOLIPs) that subsequently decay into SM states. Such particles are many cases indistinguishable 
from missing transverse energy (MET) at colliders. We propose new methods to search for these 
particles using neutrino telescopes. We study their detection prospects, assuming production either 
at the LHC or through dark matter (DM) annihilations in the Sun and the Earth. We find that 
the sensitivity for LOLIPs produced at the LHC is limited by luminosity and detection energy 
thresholds. On the other hand, in the case of DM annihilation into LOLIPs, the sensitivity of 
neutrino telescopes is promising and may extend beyond the reach of upcoming direct detection 
experiments. In the context of low scale hidden sectors weakly coupled to the SM, such indirect 
searches allow to probe couplings as small as 10~ 15 . 



I. INTRODUCTION 

Current and future experiments are beginning to probe the TeV scale. This energy scale incorporates 
the mechanism for Electroweak symmetry breaking (EWSB) and quite possibly the solution to the dark 
matter (DM) mystery. While many theoretical possibilities for TeV scale physics have been explored in 
great detail, it is still quite possible that something truly unexpected will be discovered. The reason for this 
is because no experimental evidence thus far that favors a particular model proposed so far. The possibility 
for exotic experimental signatures from beyond the SM physics has been recently explored in more detail in 
an attempt to ensure that nothing is missed at the LHC. Hidden Valleys [I], Quirks [2] and other theories 
have demonstrated the limitations of current experimental searches at colliders by exploring the possibility 
that new physics might not be triggered on, or easily reconstructed at the LHC. 

One intriguing possibility for new physics are particles which have very long lifetimes. If these particles are 
charged under the SM, they are typically known as charged massive particles (CHAMPs). These have striking 
experimental signatures, but are straightforward to search for. If the LOLIPs are neutral, experimental 
searches are much more difficult and strongly depend on their lifetime. In this paper we will explore new 
ways to search for neutral LOLIP decays that allow one to explore lifetimes that were previously beyond the 
experimental reach. 

Neutral LOLIPs are quite generic in models of physics beyond the SM. For instance, additional gauge 
sectors or new particles that couple only very weakly to the SM can easily have macroscopic lifetimes, 
r ^> O(10 _6 )s. The possibility for such lifetimes has already been discussed in many specific contexts 



2 



before. For instance, in theories of gauge mediated supersymmetry breaking (GMSB), the decay length of the 
next to lightest super-partner (NLSP) into the gravitino is macroscopic when the primordial supersymmetry 
(SUSY) breaking scale is \/F > 1000 TeV. Another scenario, that of a low scale hidden sector, has received 
considerable attention during the last year. Indeed, indirect clues from astrophysics have motivated studies 
of these "unconventional" models at the TeV scale. The recent PAMELA [3] and Fermi 4 measurements 
imply that if the experimental signals are a consequence of DM annihilations, DM must have non-standard 
properties [5], suggesting the possible existence of an additional sector, weakly coupled to the SM If this 
is the case, the low lying states in the hidden sector may be long lived and decay to the SM through some 
small coupling, e. 

Recently several ways have already been proposed to directly search for neutral LOLIP decays when the 
decay length is not much longer than the size of the detector used. These fall under the general categories of 
using high energy colliders [TUHJH], low energy e + e~ colliders [TUJ [TT] [T^] , fixed target experiments [TT1H5HT3] . 
and rare meson decays [11 . While many new ideas for searching for neutral LOLIPs have been put forth, 
the lifetimes probed are relatively short and correspondingly the coupling e is typically not smaller than 
10 -7 at most. 

An even more intriguing and less studied possibility is the case when the decay length of a neutral LOLIP 
is much greater than the size of the detector used to produce and study it. For instance if neutral LOLIPs 
are produced at the LHC and have an effective decay length significantly larger than the size of ATLAS 
or CMS, they will simply escape the detector and be recorded as MET. From this perspective one could 
not distinguish these particles from stable particles. Experimental ways to differentiate between these two 
possibilities are very limited. Two indirect methods have been suggested for probing a hidden sector that has 
particles with very long lifetimes: (i) looking at the effects on the evolution in stars, for instance supernova 
cooling [13] and (ii) measurements of certain gamma ray features in models related to DM [HI [16] . These 
methods are very model dependent and in this paper we will present more generic methods that allow one 
to search for LOLIPs with decay lengths much greater than O(km) . 

There are three ingredients which dictate the possible reach of a search for LOLIPs. The first is the 
production mechanism and its rate. This mechanism will also determine the momentum distribution of the 
LOLIPs, which may influence the capability to detect this particle. The second ingredient is the lifetime of 
the particle, which determines the probability of the LOLIPs to decay inside the detector. The physical 
quantity of interest is the effective lifetime of the LOLIPs, r e ff = jt, which depends on the typical momenta 
of the produced particle. Finally, the final states into which the LOLIP decays, will determine the types of 
detectors that are most applicable and how these particles can be searched for. 

In this paper we will focus on two mechanisms for production. The first, and most straightforward perhaps, 
is to produce the particle at a collider. If the LOLIP weakly couples to the SM but has a low mass, then in 
principle high luminosity low energy experiments can produce it [11]. Alternatively, the LOLIP could have a 
larger mass, as is the case for instance, of a bino NLSP in gauge mediation. Perhaps more interestingly, the 
long lived state might be in another sector that is only effectively accessed by producing heavier states such 
as a Z' or supersymmetric particles, which subsequently decay into this new sector. In this case luminosity 
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is still important, but higher energy is critical for producing these new states. Thus the most effective way 
for producing these particles is with the LHC. 

The second method, which turns out to be more promising, is producing LOLIPs from DM annihilation. 
For this production mechanism to be useful, one has to focus on regions where the DM density is high, 
such that a large enough flux is generated. The most interesting regions that are potentially useful are the 
Sun and Earth. Given the association with DM, there are a number of already existing constraints on this 
scenario, coming from Fermi 0], Milagro [17], and SuperKamiokande jTFj that exclude part of the parameter 
space. 

For either production mechanism, we show that using large volume neutrino telescopes is the most powerful 
method for incrcascing the reach to longer lifetimes and weaker couplings . As we will demonstrate, these 
experiments can probe hidden sectors with e as small as 1CP 15 , and lifetimes comparable to an Astronomical 
Unit, in the case of production through DM. We also demonstrate that detecting the decay of LOLIPs 
produced at the LHC is very difficult even with the improved upcoming experiments. However we suggest 
new methods that may be improved upon, and in addition put forth the idea of new detectors dedicated to 
these searches that will be investigated further in [T5] . 

The rest of the paper is structured as follows. In Section |H] we begin by discussing specific examples of 
models that contain neutral LOLIPs. We give the parametric dependence of their decay lengths, which is 
necessary to translate into specific model parameters the reaches that we present as a function of lifetime. 
In Section |III| we give the details for the neutrino telescope experiments that we will use to bound neutral 
LOLIP decays and to search for them. We include the details of how neutral LOLIP decays differ from the 



standard events looked at by these experiments and quantify this difference for certain cases. In Section IV 



we discuss potential methods for discovering neutral LOLIPs produced at the LHC. In Section [V] we discuss 
the bounds and the discovery potential for neutral LOLIP decays coming from DM annihilation in the Earth 
and Sun. In that Section, we also discuss the bounds from the Fermi and Milagro experiments in the case 
that LOLIPs arc produced from DM. Finally in Section |VI| we discuss future avenues for LOLIP searches 
and the ramifications of the bounds and reaches we calculated for individual models. 



II. MODELS 



Given that the experimental results discussed in Sections |IV| and [V] depend only on the decay length and 
the production rate, it is useful to have several benchmarks models that can accommodate the long lifetimes 
explored in this paper. In this Section we present some examples, providing the relevant formulae to relate 



these models to our general results of Sections V C We do not attempt to give a complete list of all models 
that contain LOLIPs, but the following are rather generic examples that can easily have decay lengths in 
the 1 km to 10 15 km range. 



I. Gauge Mediation. In this case the NLSP can decay to the light gravitino, G, which is typically 
assumed to be the DM particle. The lifetime is determined parametrically by the SUSY breaking 
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scale F and the mass of the NLSP, tonlspj 

/^nlspV 1 A00GeV\ 5 / y/F \ 

The lifetime can span many orders of magnitude, depending on the value of F . Unless the theory 
is complemented with a different DM candidate which is allowed to annihilate into the NLSP, long 
lifetimes in GMSB models may only be discovered using the LHC. The case where such extra DM 
candidate is introduced does not present qualitative differences to the case (Ill-d) discussed in the 
following and we will not discuss it further here. Finally, if one relaxes the requirement of a neutral 
LOLIP, another intriguing possibility present in GMSB is to have charged LOLIPs such as stau NLSPs. 
In this case very long lived staus with decay lengths of several km can also be detected in neutrino 
telescopes from production via high energy cosmic ray neutrinos |20j . 

II. R-parity Violation. R-Parity Violation in SUSY models is another typical scenario where very- 
long lifetimes may be obtained [ST]. However, as in the case of GMSB models, the R-parity violating 
MSSM does not contains long lived particles in the DM annihilation final states. It becomes therefore 
relevant for the study performed in this paper only if another DM candidate is introduced in the 
theory. Moreover, the paramctrics of the LSP lifetime strongly depends on the particular operator 
used to break R-Parity and many possibilities are present. 

III. Hidden sector with a Vector Portal. As discussed in the introduction, this theory has been 
extensively studied in the context of the cosmic ray anomalies. The hidden sector is assumed to 
communicate with the visible sector through gauge kinetic mixing with the SM hypercharge, 

evF'^ . (2) 

Here F'^ is the field strength of the gauge field(s) in the hidden sector and we have suppressed the 
adjoint index in the case of a non-Abelian group. For Abelian hidden sectors, ey is expected to be of 
order 1CP 2 — 10 since the operator above is renormalizablc and may possibly arise at the one loop 
level. On the other hand, for non-Abelian gauge groups, ey may be significantly smaller. Indeed, in 
that case, Eq. [2] arises from a non-renormalizable operator such as, 

f^V. (3) 

where M is some high scale and 4> a is in the adjoint representation of the hidden gauge group. Depend- 
ing on the ratio < <i> a > /M, ey can be as small as 10 -16 (corresponding to M at the GUT scale and 
the VEV of <& at the GeV scale). In other examples, even higher dimensional operators may produce 
miniscule couplings. We conclude that ey can span a wide range of values. 

The lowest lying state in the hidden sector may be either stable or long lived with a lifetime enhanced 
with negative powers of ey. We consider four possibilities: 

(a) Hidden photon,, jd- Id ff- If decays into the hidden sector are kinematically forbidden, the 
gauge field, 7<j, will predominantly decay into a pair of SM fermions through the kinetic mixing. 
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The branching fraction of the various channels depends on the mass of gauge fields, m 7d . The 
decay length is found to be [22] . 

cr 7W7 = (ja EM 4m 7d (l + ^ ft) * 0.8km (^) (^) , (4) 

where (3f = (1 — 4m 2 r/m 7d ) 1 / 2 is the phase space factor which we take to be 1 for the numerical 
estimate. We see that lifetimes in the range 1km < cr < 10 15 km considered here require values 
for e of order 1CT 8 1CT 15 . 

(b) Hidden Higgs, hd- hd — > //: If the hidden higgs is the lightest state in the hidden sector, it 
predominantly decays to SM fermions at 1-loop, with the approximate decay length [22 

(a d a^ M e A v m) \ / l(r 4 \ 4 /0.01 \ / m ld \ 2 / 1 GeV\ 

cr ^ * (^^< m ^ J * °- 7 x 10 km C— ; U - J (n£v) 1^7 J - (5) 

ay = .g^/47r is the gauge coupling in the hidden sector and we neglected an O(l) factor coming 
from the loop integral. Here my is the fermion mass, taken for concreteness to be the muon mass. 
As opposed to the previous case, the lifetimes of interests are obtained for relatively large values 
of ey ~ 1CP 5 -r 1CP 3 . We note, that if the higgs is not the lightest particle, but has a mass 
m ld < mh d < 2m 7tJ , it can decay through one off-shell hidden photon, in which case its lifetime is 
proportional to e v 2 but is still larger than the one of the hidden photon due to the 3-body phase 
space suppression. 

(c) Hidden Gaugino, jd- Id — ► Gl In the supersymmetric scenario, the hidden gaugino may be 
degenerate with or lighter than the hidden gauge boson. If in addition the gravitino is light, as 
in the case of GMSB, the only possible decay mode for the gaugino is into a gravitino and a SM 
photon, — > 7 + G [TSj The lifetime for this decay is 

CT _ aa tf(^y aS .&U^)'(>^'(0X. (6) 

7d ^ 7G v yiGTrF 2 ] \ e v J \ m 7i J ^lOOOTeVy v ; 

Note the large gaugino mass required to produce the lifetimes of interest. While the small mass 
limit of this model may be constrained indirectly by experiments such as Fermi |15j , large masses 
in the hidden sector are hard to probe and have not been constrained so far. The ability to probe 
such models demonstrate the strength of the indirect measurements proposed here. Finally we 
comment that this scenario can be generalized to allow other light fermions in the spectrum. 

(d) Hidden Gaugino, jd- Id — * G^d'- Another possibility is that the hidden spectrum allows the 
hidden gaugino to decay into a real hidden gauge boson. The parametrics of this decay is the 
same as in case (I) and the hidden gaugino may be very long lived. Since the long lifetime of jd is 
entirely accounted for by V~F, ev can be quite large and the subsequent 7^ decays can be prompt, 
resulting in 7^ producing final states containing SM fermions. This case presents a detection 
scenario very similar to (Ill-a), but it is achieved without the need of minuscule ey. 
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IV. Hidden sector with a Higgs portal. Another possibility for the hidden sector to communicate 
with the SM is through the so called Higgs portal, with an interaction of the form 

XH^Hh'^h' . (7) 
h' can then decay to SM fermions through the Higgs couplings. Defining the effective coupling, 

e* = A% (8) 
with v = 174 GeV the vev of the SM Higgs and its mass, the lifetime is found to be, 

' G F m 2 \ 1 /1(T 6 \ 2 ZlGeV\ 



^ = '» {^'' m '") = 30km {—) {—) ■ (9) 

where we have assumed again a decay into muons. The size of e# here is required to be small to 
allow for a sizeable lifetime. This is obtained in supersymmetric models where the interaction ^ may 
be generated in the K'ahler potential at a high scale A. In such a case, A is expected to be of order 
A ~ mf^/A 2 which can naturally be small. 

In order to detect the above particles coming from the Sun or the Earth, the DM needs to interact with 
and annihilate into the LOLIPs. In the case of hidden sector models, if the DM is not directly coupled 
to the visible sector, it interacts with SM matter and, as a consequence with nucleons, only through the 
corresponding portal. Such interactions are crucial to allow for the capture of DM in the Sun and Earth. In 
the vector portal case, the DM-nucleon cross-section is given by [BJ, 

167rZ 2 a E Madey/4„ / ev \ 2 /lGeV\ 4 
— ~ 3 x 10 cm ' 11 ' 



A 2 i 



where /i x „ is the reduced DM-nucleon mass and we have taken A = 2Z . 

Below, for models (Ill-b), (III-c) and (Ill-d), we will assume the above cross-section. While these models 
are excluded for e > 10 -6 in the case of elastic scattering of DM off nucleons [23l EI], they are not excluded 
in the inelastic case with mass splittings of order ~ 150 keV or above. On the other hand, in models (Ill-a) 
and (IV) , the relevant e is typically too small to allow for both long lifetimes and sufficiently large scattering 
cross-section (and therefore capture rates). It is possible, however, that the DM is coupled both to the 
hidden sector and visible sector with sizable annihilation fractions to both. This would be the case, for 
example, if the DM transforms under the gauge groups of both sectors. As a consequence, the DM-nucleon 
cross-section may be dominated by Z-exchange and independent of e. Below, we will assumethis possibility 
for models (Ill-a) and (IV). 



III. EXPERIMENTS AND SIGNATURES 



In this section we describe some of the experimental details needed to extract the bounds on LOLIP 
decays as well as their potential reach. As mentioned in the introduction, we focus on neutrino telescopes, 
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and in particular we consider SuperKamiokande 18J, IceCube [25 , ANTARES [26 , and the proposed km 3 
neutrino telescope in the Mediterranean KM3NeT |27j . Since we will discuss measurements from the Sun, 
we also consider bounds from the Fermi 0] and Milagro |17j experiments, which we describe in more detail 
in Section IVBl 

The probability for a particle with a decay length L = jct, to decay inside a detector of size 1 d is given 



Here D ^> d is the distance between the detector and the source. In the second equality above, we consid- 
ered the two important limits, demonstrating that if the decay length is comparable to or larger than the 
distance to the source, the probability to observe the particle increases linearly with the size of the detector. 
Conversely, if the LOLIP is short lived (compared to D) , the probability for decaying inside the detector is 
exponentially small. -Pdecay is maximized when D = L. As a consequence, LOLIPs can only be detected 
if produced at relatively nearby sources, such as the LHC, the Earth or the Sun. Looking for LOLIPs from 
DM annihilation in the Galactic center (GC) faces a suppression of at least -P<^?ay ~ 10~ 17 <i/km, which 
renders this source unusable for detection or for setting limits. Besides Pdecay the other dependence on d, D 
comes in from the usual solid angle suppression Afl — A e ff/D 2 , so that when Pdecay, is maximized the total 
suppression is approximately Vd e t/L 3 where Vdet is the detector volume. 

We focus on neutrino telescopes, since these are the largest volume detectors currently available. These 
telescopes are typical arrays of photomultiplier tubes (PMTs) in various configurations. While many different 
types of particles provide signatures in these telescopes, yet given the large volumes and coarse structure the 
events can be roughly subdivided into two categories. Muons, passing through a large part of the detector, 
leave a "track"-like signature, mostly in the form of Cherenkov radiation. On the other hand, any hadronic or 
electromagnetic particles (e.g. quarks, electrons, photons), provide more localized (typical size of O(10 m)) 
sources of light since they travel a much shorter distance than muons. With this basic categorization in mind, 
we focus on e ± or ^ ± final states as typical representatives of these two categories (we do not consider ts 
which may produce the intermediate situation of a "spot+track" ) . 

Before describing each of the experiments, a brief discussion on the case of // final states is in order and 
will apply to all the Cherenkov detectors considered below. Muon pairs originating from DM annihilations 
into LOLIPs are typically highly boosted. Given the boost factors involved, these muons will not be separated 
enough to be identified as two individual muons: for an 0(100 GeV) muon pair produced by a LOLIP with 
a mass of lGeV, the typical separation between the tracks is a few meters at most, which translates into 
O(10ns) relative time delay for the emitted Cherenkov light. Such time is comparable to typical readout 
times of PMTs. Unless specific analyses are performed, these di-muons will most of the time be recognized 



by, 




(11) 



1 The quantity d is a combination of the actual size of the detector and of the typical range traveled by the particles being 
detected, as in the case of muons. 
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as a single muon. However the Chercnkov light yield along the track of a di-muon is higher (factor of 2) than 
a single sub-TeV muon. Since the Cherenkov light yield is almost independent of the energy of the particle, 
such di-muons will be recognized as muons above their critical energy (where they start loosing energy by 
radiative processes). The critical energy in water or ice is of O(700GeV). This fact provides a good handle 
for reducing the backgrounds since the atmospheric neutrino spectra are rapidly falling with the energy. 
Moreover, by studying at the differential light yield along the track, a di-muon event will look different from 
a single muon above the critical energy: for a very energetic muon the additional light emitted on top of the 
Cherenkov emission comes from radiative processes which are most effective at the beginning of the track 
and quickly die off. On the other, hand a di-muon will yield the additional light continuously along the 
track, providing another possible handle for isolating the signal. It would be very interesting to study this 
effect further and the implications it may have for energy measurements in neutrino telescopes |28j . 

Turning to the specific experiments, SuperKamiokande is a large volume of water surrounded by PMTs. 
We can use its data from upgoing muons to constrain LOLIPs from the Earth and the Sun. These LOLIPs 
may either decay into muons in SuperK, or decay earlier, producing neutrinos via subsequent muon decays. 
These interact with the material and produce muons that will be detected. It is important to note that 
looking at the combined total flux provides only a rough bound because, as explained above, the events 
that come from LOLIPs which decay inside the detector are quite different from standard charged current 
neutrino scattering. 

The ANTARES experiment is a more modern neutrino telescope and is a test project for a future km 3 
array in the Mediterranean, KM3NeT. While this experiment is not the largest neutrino telescope, it is 
located approximately 300 km from the LHC which makes it the most compelling detector when discussing 
LOLIP production at the LHC in Section |IV| The experiment is comprised of 12 strings at a depth of 
approximately 2500 m underwater in the Mediterranean sea. Each string has 25 storeys separated vertically 
by 14.5 m and the strings are spread out horizontally from each other by approximately 70 m. Each storey 
contains an array of three large PMTs. The effective area for this experiment is ~ 0.1km 2 for very high 
energies and drops rapidly towards the 100 GeV range where the effective area is of order of 0.02 km 2 . For 
a LOLIP that decays into fi + fi~, we have to take into account that the muon pair has a very small opening 
angle and there therefore an increased light yield which in turn affects the effective area, A e ff, compared to a 
single muon event. We have done a Monte Carlo simulation of this effect taking into account the ANTARES 
geometry and trigger table (as well as other effects like the PMT response as in or the water attenuation of 
light) and we nominally find that the effective area of ANTARES for a O(100 GeV) horizontal di-muon is 
roughly 3 times larger than the case of a single muon at trigger level. 

When discussing the potential reach for neutral LOLIPs, we will also examine the planned KM3NeT 
telescope. This turns out to be more interesting for potential LHC LOLIP production if KM3NeT is built at 
the site of ANTARES. KM3NeT is a proposed km 3 neutrino telescope to be built at one of several locations 
in the Mediterranean sea. To determine its possible effectiveness, we used the proposed ANTARES-like 
design for it. This consists of a homogeneous array of 225 strings separated by 95 m from one another, with 
37 storeys each spaced 15.5 m apart. For simplicity we use the original optical module design of ANTARES 
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and the same triggers to MonteCarlo the detector response. We find that for muon energies of ~ 100 GeV 
the effective area for KM3NeT is about a factor of 30 larger for LOLIP decays than ANTARES at the same 
energy. The KM3NeT effective area also sharply falls below 100 GeV based on our simulations, the main 
reason being the large spacing between the strings. 

The largest available neutrino detector currently is the IceCube detector. IceCube [25] is a km 3 detector 
which consists of 80 strings at depths between 1400 m and 2400 m under the Antarctic ice. The strings 
are distributed over an area of approximately 1 km 2 and each contains 60 optical modules spaced 17 m 
apart. For this experiment we use the effective areas given in We also investigate the inclusion of the 
DeepCore extension (a more densely instrumented inner volume designed to reduce the energy threshold 
of IceCube) , by correcting these effective areas according to the preliminary information of [30] . This does 
not take into account the effects associated with the special type of events that LOLIPs produce with two 
highly collimated charged particles. The reach we give in Section [V] will thus be conservative. It would be 
interesting to determine exactly how sensitive neutrino telescopes are to these types of events. 



IV. LHC PRODUCTION 

When LOLIPs arc produced at the LHC, some fraction of them will decay within the volume of the ATLAS 
or CMS detectors. Given a canonical size for these detectors of approximately d ~ 10 m, and assuming a 
decay length L, the amount of particles that decay within the detector is given by, 



^decayed — ^produced 



(l-e- d ' L ), (12) 



which, as in Eq. (Ill gives 



^decayed ~ ^produced j^- (1^) 

This also holds if only part of the detector can be used for detecting the LOLIP decays, as in the more 
realistic size of lm. Taking a typical strong production cross-section of order O(100 pb) at the LHC for 
heavy particles, one obtains N proc i U ced ~ 10 6 events with an integrated luminosity of 10fb _1 . Consequently, 
for a decay length of order a few hundred meters, thousands of LOLIPs may decay inside the LHC detectors. 
There has been a great deal of work on the possibility of searching for highly displaced vertices at both 
ATLAS and CMS. However, given the many backgrounds, it is not clear how many events one would need 
for the LHC detectors to be useful in studying neutral LOLIP decays. Moreover, for lifetimes comparable 
or greater than O(100 km), the number of events at the LHC becomes negligible and the LOLIPs will only 
show up as MET. 

In order to differentiate between stable particles and LOLIPs external detectors are needed. The most 
straightforward idea would be to use other detectors to conduct a long baseline experiment similar to MI- 
NOS [31 . The nearest large volume detector to the LHC is the ANTARES neutrino telescope in the 
Mediterranean sea. A rough estimate of the detection rates at ANTARES is given by the number of par- 
ticles decaying inside the ANTARES detector. If we assume the effective lifetime of these particles is of 
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0(300 km), and assuming they are sufficiently energetic such that A e ff ~ 0.1km 2 , one has 




'produced ~ 30047r300 2 




(14) 



This means that for cross-sections below O(100 nb) ANTARES will not be useful for detecting such LOLIP 
decays. If the KM3NeT experiment was located at the same site as ANTARES, it would in principle have 
a better chance to discover LOLIP decays. To illustrate the utility of KM3NeT we choose a benchmark Z' 
candidate with a 500 GeV mass that decays into two LOLIPs which subsequently decay into a pair of muons. 
The Z' model point we chose has a production cross-section of 12 pb. If we choose a lifetime cr = 10 5 cm to 
minimize the suppression of the decay probability, we find that for the optimal orientation of the KM3NeT, 
1.5 x 10~ 2 events will be observed per year. This should be compared with 6 x 10~ 4 events will be observed 
per year for ANTARES. The event rate for this benchmark is clearly too low. Other models such as GMSB 
where one could use the entire SUSY cross-section would be more promising, however with the current LHC 
designed luminosity, this is doubtful as well since KM3NeT will only be sensitive to O(nb) event rates. It 
is possible that Super-LHC (for a review see 32J) may be sufficient to allow for a detection of LOLIPs at 
KM3NeT. We postpone this question to future work. 

The low number of events can be understood by noting the tension between the need to raise the collider 
energy in order to increase the effective area and the need to lower the energy in order to increase the 
production rate. It would therefore be advantageous to lower the threshold in KM3NeT as much as possible. 
Despite the low event rates, it is expected that the background in these neutrino telescopes can be removed 
entirely with the use of angular information and timing synchronization with the LHC detectors. Thus 
observing any O(l) number of events could be a clear signal of LOLIP decay. 

Given the small rates associated with the nearest neutrino telescopes one is forced to examine other 
alternatives for studying LOLIP decays produced at the LHC. The suppression due to the solid angle clearly 
indicates that new detectors that are closer to the LHC are needed. Additionally one could attempt to use 
the existing LHC detectors as long baseline detectors for each other. Events at ATLAS could be recorded 
at CMS or vice versa. These possibilities require detector design studies, and a more detailed analysis of 
timing between the experiments and the triggers. We leave both possibilities for future work [19j . 



In this Section we examine the prospects of observing LOLIPs which are produced from DM annihilation. 
The most promising sources to look for are regions where the DM density is enhanced, such as the Sun, 
Earth, sub-halos and Galactic Center (GC). Out of those, only the Sun and Earth are suitable for detections 
in neutrino telescopes since the probability of LOLIPs to decay inside the detector is too small for more 
distant sources. Below we consider the constraints from SuperK, Fermi and Milagro on LOLIPs originating 
from the Sun and Earth, and study the reach of IceCube for discovering these particles. We begin by shortly 
reviewing the formalism and relevant constraints before presenting our results. 



V. 



LOLIPS FROM THE SUN AND EARTH 
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A. Production Rate 

The production rate of LOLIPs in the Sun and Earth depends on the capture rate, which can vary by 
many orders of magnitude depending on the nature of the DM. This topic has been studied extensively in 
the literature j33j[34] and here we shortly review the necessary ingredients. 

To determine the rate one solves the evolution equation for the DM number density, 

N = C — C A N 2 . (15) 

Here C is the capture rate and Ca = {, av )l^eSf ^ s the thermally averaged annihilation cross section, 
per effective volume, with V eff = 1.8 x 10 26 (m DM /TeV)" 3 / 2 cm 3 for the Sun [31 |35] and V eff = 
5.7 x 10 22 (mDM/TeV) _3 / 2 cm 3 for the Earth (33j [34]. Solving the above equation, the DM annihilation 
rate in the Sun is found to be, 

r= \c A N 2 = ^Ct<mh 2 {y/CCA~t), (16) 

where t = 4.5 x 10 9 years is the age of the solar system. For sufficiently large capture rate and DM 
annihilation rate, the DM number density saturates and the annihilation rate can be well approximated 
by T = (7/2. This approximation holds for C > C§ = 3 x 10 17 (m DM /TeV)~ 3 / 2 sec" 1 for the Sun and 
C > Cm = 10 14 (rnDM/TeV)~ 3 / 2 sec -1 for the Earth, where we assumed the usual WIMP annihilation cross- 
section, (civ) = 3 x 10~ 26 cm 3 sec -1 . This approximation turns out to be rather good in the Sun but not in 
the Earth where equilibrium is typically not reached. 

The above solution depends both on the DM capture and annihilation rate. The former is model dependent 
and in particular depends on whether the DM interacts clastically with the nucleus as assumed in many direct 
detection searches. If there are more complicated interactions with nucleons, as in Inelastic DM [36 or in 
other recent models |37j . the capture rate can be significantly modified. In order to demonstrate the effects of 
a modified capture rate, we consider elastic scattering and inelastic scattering. For the elastic case see [531134"] 
while the inelastic case has been recently studied in [35J I n the inelastic scenario, the DM states are 
split with a mass difference 8, and scattering off nuclei excites the light state into the heavier one. While 
the net effect of inelastic scattering is to lower the capture rate, its virtue is in relaxing the constraints on 
the DM-nuclcon cross-section arising from direct detection, by as much as three orders of magnitude. This 
allows for a significant enhancement in the annihilation rate from the Sun. On the other hand the rate from 
the Earth shuts down already for very low mass splittings, since elastic scattering always dominates there. 

We have repeated the analysis of [38, 39 and have reproduced both the elastic and inelastic capture rates. 
As a reference, it is convenient to express the capture rate as a function of both the DM mass and splittings 
in a closed form. We find, 

r ffx , /!TeV\ 3/2 / p DM t \ / 250kms^ \ 3 / g* \ 
\ m DM / \0.3GeVcm A J \ Vdisp / \ 10 cm / 

with i = ©,© refers to the Sun and the Earth cases, pdm is the DM density in the solar system, Wdisp is 
the DM velocity dispersion and <7gj is the spin-independent DM-proton cross-section which we assume (as is 
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typically the case) to dominate over the spin-dependent cross-section. The constants have values 

rc ~ 2.2 x 10 20 s _1 « e ~ 6.6 x 10 10 s _1 . (18) 

The functions fi are approximate functions fitted to the inelastic results, and are given by, 

tanh(l-V(120keV)) + l 
/oOWdm) = tanh(1) + i • (19) 

fm(S,m BM ) = 8(6-0 keV) (20) 

The above fit is correct at the 20% level for ttidm ^ 200 GeV and to about 50% at m^M = 100 GeV. We 
note that inherent uncertainties arising from the astrophysical uncertainties are significantly larger, thereby 
rendering the above approximation sufficient. 

B. Existing Constraints 

While a great deal of the parameter space for neutral LOLIPs is unexplored, some constraints already 
exists. Collider constraints from low energy experiments have been covered in many different ways in [7J 
HJ [lOj [TTJ Q21 [13l [14]. However, if DM particles can annihilate into LOLIPs there are additional indirect 
constraints on some part of the parameter space. These constraints roughly fall into three categories: existing 
neutrino telescope bounds, photon and electron fluxes from the Sun, and DM direct detection experiments. 
If the LOLIPs subsequently decay into muons, then SuperK bounds from up-going muon searches from both 



the Sun and the Earth are relevant [IH1 HI] . As explained in Section III these searches are just a rough 
bound on the LOLIP scenario because utilize on neutrino induced muon events in SuperK. 

As long as the LOLIP effective lifetime is long enough such that they escape the Sun, it will produce 
decay products between the Sun and the Earth. These decays may give rise to signals which can be explored 
using experiments other than neutrino telescopes. For instance if the LOLIP decays to charged particles 
these particles will radiate photons on their way from the Sun to the Earth that can be detected by gamma 
ray searches such as Fermi or Milagro. If the LOLIPs decay products ultimately produce e , additional 
bounds coming from the Fermi measurement of the total (e + + e~) flux apply. Specifically for the Fermi 
7 measurements we consider the differential 7 ray spectrum for the Sun, as was recently presented by the 
Fermi Collaboration jJD], extending up to 10 GeV and fitted by a power law with a spectral index of -2.25. 
We require that the LOLIP contribution should not exceed twice the size of the error bars in any given bin. 
The same procedure is applied for the case of (e + + e~) Fermi data. On the other hand, the Milagro 7 data 
extends to much higher energies than those of Fermi. In order to compute the bounds we use the publicly 
available effective area and the exposure time to extract the total number of events expected in Milagro. 
We then extract the 90% CL intervals according to the procedure described in [17]. Before passing to the 
Direct Detection experiments, we further note that even though these bounds are directly applicable to the 
Sun, they also indirectly bound the production of LOLIPs in the Earth because the capture and annihilation 
rates in the two astrophysical objects are related. 
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FIG. 1: The IceCube reach for detecting LOLIPs produced in the Sun and decaying into muons, as a function of the 
annihilation rate and decay length. On the x-axis, -Rhid — <r(DMDM — > LOLIPs) /ototai and BR^ is the branching 
fraction of LOLIPs into muons. 100 GeV and 2 TeV DM masses are shown, assuming DM interacts inelastically with 
mass splitting of 200 keV and decays into 500 MeV LOLIPs. The shaded regions are excluded by SuperK (yellow), 
Fermi photon measurements (red) and Fermi (e + + e~) measurement (orange). The blue and green solid lines show 
the one and five year sensitivity of IceCube to measuring LOLIPs. To demonstrate the sensitivity to the special 
di-muon events from LOLIP decays inside the detector, we show in dashed lines the one and five year discovery reach 
for such events only. For reference two theory lines corresponding to models (Ill-a) and (Ill-b) described in Section |Tl| 
are shown. Finally the dashed gray line is the current CDMS bound on the annihilation rate. This constraint is 
irrelevant in the 100 GeV case. 

The last type of bound from DM production of LOLIPs is from DM direct detection experiments. This 
bound comes from the fact that any capture of DM in the Sun or Earth than requres a non- vanishing DM- 
nucleon interaction. Using the relation between the annihilation and capture rates we can then translate 
a point in the parameter space of lifetime and annihilation rate into a specific value for ct x „. Requiring 
that this interaction cross-section does not violate the constraints from the XENON10 [24] and CDMS [23] 
experiments, one can bound the parameter space. It should be noted that this translation is model dependent, 
e.g. it depends whether the scattering is elastic or inelastic. Recently there have been many other models 
proposed that can also change the nature of DM direction detection [37] : We do not analyze the bounds for 
all these cases, but warn the reader of the many caveats in the existing direct detection constraints. 
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FIG. 2: Similar to Figure [JJ for the case of annihilation in the Earth. The inelastic splitting is taken to be zero since 
the capture rate is highly suppressed otherwise. The red shaded region shows the constraints from the Sun, assuming 
that a non-vanishing annihilation rate at the Earth implies a related annihilation rate in the Sun. 



C. Results 

1. Sun 

Let us now present the results for the discovery reach of LOLIPs from the Sun. In Fig. [l] we show the 
constraints which arise from SuperK, Fermi and CDMS together with the detection reach in IceCube for 
one and five years, in the case of a LOLIP decaying into muons. These plots and the ones to follow show 
the sensitivity of IceCube in the decay length - annihilation rate plane. The measured rate is suppressed 
with respect to the annihilation rate, by Rud — <r(DMDM — > LOLIPs)/cr t otai and by the specific branching 
fraction of LOLIPs studied. Before proceeding, a few words on our treatment of the backgrounds are in 
order. 

The IceCube curves represent the 95% CL limits. We use the atmospheric neutrino flux [IT] to estimate 
the background. For the case of LOLIPs decaying in proximity of the detector, we considered the expected 
atmospheric flux above 500 GeV, independently of the LOLIP energy, because to the different Cherenkov 
light yield of a di-muon pair. On the other hand, for the events initiated by neutrinos from earlier muon 
decays, we consider the full atmospheric muon background above the detection threshold. We study the flux 
coming from a cone of 0(3°) around the source of interest. In the case of SuperK, we directly use the bounds 
provided for standard DM i/„ searches [T5] and also consider the showcring-muon analysis [35] to directly 
bound LOLIPs decays (including also the showering-/^ fraction from events as outlined in [25] S3)- 

In both plots of Figure [l] we have taken the LOLIP mass mLOLiP = 500 MeV. The main features of the 
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figures are easy to understand: The detection rate is maximized when the LOLIP decays at or just in front 
of the detector. Given the LOLIP boost of order 7 ~ wdm/^lolip, one finds a linear decrease in the reach 
for lifetimes above the critical value, ct'' ~ iulolip ' 1AU/jt»dm = 0.5GeV AU/moM; where AU= 1.5 x 10 8 



km. This linear decrease is in accordance to Eq. (11). Below cr c , we find that the sensitivity is roughly 
constant as long as the LOLIP escapes the Sun. This can be understood as follows. LOLIPs that decay 
to muons before reaching the Earth produce neutrinos. The flux of these neutrinos is dominated by the 
decays occurring very close to the Sun's surface and the probability of detecting them is suppressed by the 
probability of the neutrinos to interact with the detector, 

P„-- mt ~ (R M + d)x P N A x a vN . (21) 

Here d is the length of the detector, is the muon range, p is the density of matter in the detector, (ice 
or water), Na is the Avogadro number and o v n is the Neutrino- nucleon cross-section. Rather amusingly, 
Pv-int ~ -Pdccay for mrjM — 100 GeV, which explains the flat region in the plot as ct drops below cr c . The 
flat region stops at the value of ct corresponding to the Sun radius. The products of LOLIPs that do not 
escape the Sun are not detected, as in this case the muons stop inside the it giving rise to neutrinos of very 
low energies. In the plots, the dashed lines show the sole contribution of muons from LOLIPs decaying nearby 
the detector, while the solid lines include the additional neutrinos produced before reaching the Earth. As 
discussed in Sec. |III| the tracks of di-muons produced inside the detector have characteristic features which 
should allow one to differentiate them from regular single muon events. It is therefore conceivable that one 
can further reduce the backgrounds improving these results. Moreover, by measuring both the total flux and 
the di-muon fraction, a rather good estimation of ct may be feasible. We postpone this analysis to future 
work. 

Figure[l]also shows two calibrated "theory" lines, corresponding to models (Ill-a) and (Ill-b). As discussed 
in Section [TT| Model (Ill-a) describes a hidden sector with a long lived hidden gauge field which mixes with 
the SM through gauge kinetic mixing. We see that IceCube can probe extremely small values of ey- Unless 
the DM interacts directly with the SM, the capture rate and hence the annihilation rate from the Sun 
vanishes for such low values. Model (Ill-a) therefore assumes that the DM is also coupled to the SM (see 
Section |n| and we assume a DM-nucleon cross-section corresponding to an annihilation rate of 10 21 sec _1 , 
below the CDMS bound. The same situation also applies to model (IV), hence this line can also be used to 
read off the numbers for this model as well, upon a trivial rescaling of en — 6 • 10 2 ey- For the other theory 
line as in model (Ill-b), the lowest state is a Higgs which can decay to SM fermions at one loop. Here, the 
lifetime is proportional to £y 4 and hence IceCube probes rather large values of ey. 

In Figure. [3^ we plot the case for jtidm = 2 TeV where LOLIPs decay predominantly into electrons. This 
would occur when the LOLIP is very light, below twice the muon mass. We therefore take tololip = 25MeV. 
As discussed in the previous Section, the detection of electrons is different since they interact with the 
detector promptly to produce a localized signal rather than a long track. Therefore the analysis should also 
be indicative of the reach for other final states producing electromagnetic or hadronic showers. Unlike the 
muon case, we do not consider the SuperK bounds since the photons constraints from Fermi are comparatively 
much stronger than in the muon case. This is partly due to the fact that the background from atmospheric 
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FIG. 3: Left: The reach for inelastic DM annihilating in the Sun into LOLIPs with a mass of 25 MeV, which 
subsequently decays into e + e~ . For reference two theory lines corresponding to models (Ill-a) and (Ill-b) described 
in Section |ll] are shown. In addition to the constraints shown in Figure. [I] we plot the constraint on the photon flux 
from Milagro which is insignificant in the muon case. Right: The reach for annihilation of DM in the Sun into a 
LOLIP with a mass of 50 GeV, subsequently decaying via hard photon as in Model (III-c). The axes in both planes 
are as in Figure [T] 

electron neutrinos is larger by roughly two orders of magnitude [43] than in the case of muons because of 
poorer directionality reconstruction for showers compared to /i tracks. Unless a more detailed analysis is 
performed and better ways of background rejection are found, the bounds from Fermi gamma ray searches 
cover almost the full reach of IceCube and no discovery can be made in this case. In Fig. [SJa we also plot the 
IceCube reach for the case where LOLIPs decay into hard photons, and include a theory line for the model 
(III-c) described in Section |ll] As one can see, Milagro gamma rays constraints cover most of the IceCube 
reach. 

2. Earth 

Let us now turn to the Earth. Due to the different scales involved, the search for LOLIPs arising from 
DM annihilations at the center of the Earth enables one to probe different lifetimes, of order cr ~ R^/m^M- 
As opposed to the Sun, LOLIPs that decay before reaching the detector stop instantly and therefore if any 
flux is measured, it must be directly related to LOLIPs decaying nearby IceCube. Furthermore, as discussed 
earlier, the capture rate drops quickly as the mass splittings between DM states is increased. Consequently, 
we concentrate on the elastic scattering region. In Figs. [2] and [4] we present similar plots to those of the 
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FIG. 4: Left: As in the left panel of Figure[3j only computed for annihilation in the Earth, no inelastic scattering, and 
DM mass of 2 TeV. Right: Similar to the right panel of Figure [3] The red region in both plots show the constraint 
from the Sun, assuming that a non-vanishing annihilation rate at the Earth implies a related annihilation rate in the 
Sun. 

Sun. We take 6 = OkeV and plot the reach for IceCube as well as constraints from CDMS and SuperK 
(for the muon case). We also plot the strongest constraint from the Sun, which corresponds to the Fermi 
and Milagro bounds on the gamma ray flux. This assumes that the capture in the Earth would occur in 
conjunction with capture in the Sun. As before, in the case of e and 7 final states we do not include the 
SuperK bound for decays into electrons. This is because of the much larger backgrounds compared to muon 
searches, that weakens the bounds. It is apparent that the allowed region that can be probed by IceCube is 
significantly smaller for the Earth case. Nevertheless, such measurements allow us to cover a different region 
of parameter space. 



3. Comparison to Direct Detection Bounds 

It is interesting to understand whether the region in reach of IceCube can be ruled out by direct detection 
experiments in the near future. The answer is presented in Fig. [5j where we plot the maximal sensitivity of 
IceCube to the DM-nucleon cross-section as a function of the DM mass. To extract this plot we pick, for a 
given DM mass, the optimal lifetime that would allow for the lowest annihilation rate. We translate it to 



the nuclear cross-section, assuming a particular DM annihilation cross-section, as outlined in Eqs. (17l,(18l. 
The curves are drawn for the final states both for the case of the Sun and the Earth. For concrctcncss we 
fixed the LOLIP mass to 500 MeV. Since we scan over the LOLIP lifetime to achieve maximal sensitivity, 
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FIG. 5: In both the left and right panels of this figure we plot the reach for IceCube in the standard m x vs a xn 
plane used for plotting constraints from direct detection experiments. The procedure for translating from T and ct 
to m x and a xn is described in the text. Additionally we plot the current bounds from CDMS and XENON10 and 
the projected sensitivities for upcoming DM direct detection experiments. 

the effect on Fig. [5] of changing mLOLiP is only subleading. 

The solid blue(green) line corresponds to the one(five) year reach for IceCube in the case of a WIMP 
annihilation cross-section, a = 3 x 10 _26 cm 3 s _1 . The dashed lines correspond to a = 3 x 10~ 23 cm 3 s 
which is the cross-section required to explain the recent cosmic-ray anomalies [5] . Nucleon scattering cross- 
sections of 10~ 46 cm 2 can be achieved in 5 years in the case of the Sun irrespectively of the annihilation 
cross-section, while the Earth can reach that level only if DM annihilation is enhanced [44]. We compare the 
results with the given and future constraints from direct detection. As one can see, for the case of the Sun, 
IceCube comparatively performs better at higher energies. We stress that the comparison is made for the 
elastic, 5 = keV, case. For inelastic scattering, constraints from direct detection are significantly weaker 
and the corresponding IceCube reach (modestly weaker as well) can be easily inferred from the approximate 
formulas we have given. Therefore, we conclude that IceCube will have a reach better or comparable to that 
of direct detection experiments and will provide an additional avenue to test non-standard DM scenario. 

4- Constraints on Models 

Finally we conclude this Section by presenting the implications of the results described here for the models 
of Section |nj This is illustrated in Figure [6] for the models (III-a,b,d) and (IV). The experimental reach and 
the current constraints in terms of the fundamental parameters of each model are shown. On the top-left 
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FIG. 6: The current bounds and the IceCube reach in the parameter space of four benchmark models presented in 
Section [Tl] are shown based on the results summarized in this section. Upper left: model (Ill-a). Upper right: model 
(Ill-d). Lower left: model (Ill-b). Lower right: model (IV). The red band correspond to the photon exclusion region 
extracted from Fermi and Milagro data. 



panel we plot the Sun reach for a 2 TeV DM annihilating into a vector portal described in model (Ill-a), 
assuming mass splittings of 200 keV and ey-independent annihilation rate of 10 21 s . On the top-right panel 
we show a similar plot for the vector portal model, (Ill-d), where the LOLIP is the gaugino with a lifetime 
controlled by the SUSY-breaking scale. The left-bottom panel is for the 2 TeV DM with the light higgs in 
the vector portal, and finally the bottom-right is the higgs portal, model IV. In all plots the excluded part 
of the parameter space due to the Fermi photon measurement is shown by the red shaded region. 
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By changing the mass of the particle decaying back to the SM (which, except for model (Ill-d) corresponds 
to the LOLIP mass) the branching fractions into final state particles are changing. In particular in the case 
of models (Ill-a) and (Ill-d) the final states are selected by the requirement that the mediator couples to 
the electric charge, while in the other two cases by the requirement that the partial widths are proportional 
to mj. For mediators lighter than 1 GeV this results in a combination of e + e _ , ^t + /i~ and tt + it~ , 7r°7r° final 
states. We estimate the admixtures by the measurement of the e + e~ — » hadrons cross-section for the case of 
the electric charge, while in the other models we estimate the contribution of pions using chiral perturbation 
theory. All the kinks and drops shown in these plots can be readily understood in terms of the relative 
importance of these channels. It is notable that we do not extend the plots for masses above 1 GeV (where 
low energy direct probes lose their sensitivity). The reason for this is that while allowed regions may still 
exist up to masses of the order of few GeV (photon constraints definitely becomes very stringent somewhere 
before 8 GeV [5]), QCD uncertainties do not allow for a satisfactory theoretical or experimental handles to 
determine the possible SM final states and their relative branching fractions. 

VI. CONCLUSIONS 

In this paper we have studied the prospects of detecting neutral long lived particles (LOLIPs) with decay 
lengths ranging between 1 km and 10 15 km. Such long lived states may arise in several scenarios such as 
GMSB models and hidden sectors weakly coupled to the SM. Previous studies focused on significantly 
shorter lifetimes, keeping many theoretical models out of the experimental reach. We stress that the added 
value of the approach discussed in this paper, is that it not only improves on the sensitivity to low scale 
weakly coupled sectors, but also it allows for detection of hidden sectors with mass scales significantly larger 
than 1 GeV, such as the model (Ill-d) of Section [TT] 

To obtain sensitivity to such long lifetimes, we studied neutrino telescopes employing their large scale 
detectors. For the production mechanism, we considered the LHC and DM annihilations into LOLIPs in 
the Sun or the Earth. Attempting to directly probe LOLIPs produced at the LHC and detected at Antares 
or KM3NeT proved not feasible due to the low rates. Future improvements on the luminosity at the LHC 
and/or on the energy threshold at KM3NeT may render this possibility viable. Conversely, designing small 
scale detectors closer to the LHC may allow one detect LOLIPs which show up as MET at the LHC. We 
postpone the study of this possibility to future work [TIT]. On the other hand, indirect detection of LOLIPs 
from the Earth and Sun seems promising, depending on the final states into which the LOLIPs decay. 
In particular, if LOLIPs predominantly decays into di-muons, a significant region of the parameter space 
previously unexplored can be probed by IceCube. Conversely, the prospects of detecting final states which 
shower inside the detector, such as electrons or photons, strongly depend on the capabilities of reducing the 
atmospheric neutrino background. 

Finally, in relation to DM annihilating into LOLIPs, we have compared the capabilities of the neutrino 
telescope techniques, to those of direct detection. We found that for models which accommodate long lived 
states, neutrino telescopes can detect annihilation fluxes from the Sun corresponding to nucleon cross-sections 
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beyond the reach of upcoming direct detection experiments. 

Note Added: While this work was in writing, related work appeared [IS]. 
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